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Abstract 

Single-particle spectra of A and S hypernuclei are calculated 
within a relativistic mean-field theory. The hyperon couplings used 
are compatible with the A binding in saturated nuclear matter, 
neutron-star masses and experimental data on A levels in hypernu- 
clei. Special attention is devoted to the spin-orbit potential for the 
hyperons and the influence of the /9-meson field (isospin dependent 
interaction) . 

PACS number: 21.80.+a 
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Relativistic Mean-Field Calculations of A and S Hypernuclei 



N.K. Glendenning 
D. Von-Eiff, M. Haft, H. Lenske, and M.K. Weigel 

Nowadays the hypernuclear physics are of great interest in many branches 
of physics. Of particular importance is the understanding of strange parti- 
cles in baryonic matter, since many questions in heavy-ion physics and as- 
trophysics are related to the effect of strangeness in matter. Experimentally 
and theoretically such problems have been mainly studied for the A hyperon, 
because it has zero isospin and charge. For A particles theoretical approaches 
for the spectroscopy range from nonrelativistic models to the relativistic 
Hartree approximation (RHA) 0]. For the RHA one uses a Lagrangian 
with effective A couplings to the a- and w-meson fields. In Ref. the A 
coupling constants (i.e. their relative strength to the corresponding nucleon 
couplings Xfj = Qha/ga and = gKuilgui) have been fitted to the experi- 
mental A hypernuclei spectra. However, treating and x^^ as independent 
parameters leads to a highly uncertain determination (correlation errors up 
to ±65% in Ref. g). 

On the other hand, the contribution of the hyperons strongly influences 
the mass of neutron-stars. In a recent publication 0| Glendenning and 
Moszkowski related the scalar and vector couplings of the A hyperon to its 
empirical binding in saturated nuclear matter 0] and thereby obtained com- 
patibility of this binding energy with maximum neutron-star masses. In fact, 
the large correlation found in the least squares fit mentioned above Q reflects 
this relation of x^ and x^^ to the A binding in nuclear matter. In summary, 
one finds that (1) neutron-star masses, (2) the A binding in saturated nuclear 
matter, and (3) A levels in hypernuclei are mutually compatible and rather 
narrowly constraining the A couplings. 

Concerning S hypernuclei, up to now the experimental situation is not 
satisfactory. S hypernuclear production has been investigated at CERN, and 
later at Brookhaven and KEK, but the statistical accuracy of the available 
data is not very good, because of the strong conversion the S undergoes in the 
nucleus {T,N —>■ AN). The controversial evidence for narrow S states (F < 
5 — 10 MeV) is reviewed in Ref. Therefore, in theoretical investigations 
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including S hyperons one usually assumes an universal hyperon coupling; 
i.e. all hyperons in the lowest octet have the same coupling as the A 0, 
1^. The prospects for significant advances in high resolution hypernuclear 
spectroscopy at CEBAF or at future facilities such as the proposed PILAC 
and KAON are discussed in Ref. [0. 

It is the aim of this contribution to analyze A hypernuclei under consid- 
eration of the constraints (l)-(3) mentioned above and to extend such an 
investigation to S hypernuclei. 

For the nucleonic sector we use the common nuclear field theory La- 
grangian including the nucleon couplings to the a-, cu-, and p- meson fields 
[§] plus phenomenological cr-selfinteractions For the three charge states 
of the S hyperon we write the following Lagrangian WU, [TTI]: 



s 

- S.,f^7.e;, + ^,A-(r3),,)s,„ (1) 



where Ejj and 0^- are the traceless 2x2 matrices 



and 

~ [ V2p^^ -P'o I • 

The sum on S in the first line of Eq. (|l]) is over the charge states S~, E°, 
and S+. The Euler-Lagrange equations then yield the Dirac equations for 
the S hyperons: 

{i^^d^ - M* - gj^^^^uj^ - h^gj^plf^P^ - /sscTm^'^) = 0, (4) 

M* = Me - gsa(r, (5) 

where /3s denotes the third isospin component; i.e. /3s = — 1,0,+1 for S~, 
E°, and E"*", respectively. This means, that Eq. (Q) already considers the 
fact, that within a RHA only the charge neutral component of the p-meson 
field, Po, yields a nonzero ground-state expectation value. 



2 



The A hyperon has isospin and charge zero and therefore cannot couple 
to the p-meson and electromagnetic fields. Hence, under the consideration of 
the universal hyperon coupling and the replacement of Ms by Ma in Eq. (j^), 
the A Dirac equation equals the S*^ Dirac equation. 

For calculating the hypernuclear spectra we made use of a technique sim- 
ilar to the so-called expectation-value method, which was successfully used 
within nonrelativistic and relativistic nuclear physics to incorporate shell ef- 
fects into semiclassical densities and energies (see Refs. [12-14]): The Dirac- 
Hartree equations for the hyperons are solved only once with the meson fields 
of the corresponding nucleonic system, which are self-consistently determined 
within a relativistic Thomas- Fermi (RTF) approximation, as an input. 

To check the validity of this approximation we recalculated the A single- 
particle spectra for the hypernuclei A'Ca and "j^^Ph with the parameters of 
Ref. 1^ (parameter set I in Table I). The results for various A levels are 
displayed in Table II where our results are denoted by H*. Compared with 
the fully self-consistent RHA results the H* approximation systematically 
underestimates the A bindings which may be attributed to a surface energy 
that is somewhat too large within the RTF approach [T^. But as expected, 
the agreement is better for the larger mass number A and the deeper lying 
levels because in both cases the Thomas- Fermi assumption of locally constant 
fields is more valid. In conclusion, our results show a rather good agreement 
with those of Ref. 0, which gives confidence in the described scheme. 

In the next step we calculated several A and S hypernuclei using a set 
of coupling constants from Ref. ||^, which considers the constraints (l)-(3) 
mentioned above and in addition has been successfully used in the description 
of nuclear matter properties (parameter set II in Table I). In Fig.l we show 
the contributions of the meson and electromagnetic fields to the hyperon 
self-energy for the nuclei ^^Si, ^°Ca, ^''Zr, and ^°^Pb. The nonrelativistic 
reduction of the hyperon potential, (for A and S° entirely, for mainly) 
given by the difference gnuj^^ — gHa<^, if = A, S, is also displayed. It is 
worth noting, that the small potential depths of ~ 30 MeV go along with 
a relatively smooth radial dependence (compared with nucleonic potentials), 
thereby additionally supporting the feasibility of the Thomas-Fermi meson 
fields. A similar behaviour was found for O within the RHA calculations 
of Ref. 0. In Figs. 2 and 3 we show the single-particle spectra of protons, 
neutrons. A, E*^, and S"*" hyperons for the nuclei ^°Ca and ^°*Pb. Because 
of the smaller couplings the hyperon levels are considerably less bound than 
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the corresponding nucleon levels. Looking at the A and S*^ single-particle 
energies, the larger yields a smaller repulsive effect of the kinetic energy 
resulting in systematically stronger bindings for the S". 

Dealing with hypernuclear states and their structure, one of the most in- 
teresting questions concerns the spin-orbit potential for the hyperons 0. It 
is one of the great advantages of a relativistic treatment that the spin-orbit 
interaction is automatically included in the single-particle Dirac equation, 
and can be identified by means of a Foldy-Wouthuysen reduction. For exam- 
ple, looking at the charge neutral A and S° hyperons, both of which are not 
coupled to the p-meson field, the ratio of the spin-orbit splitting (Thomas 
terms) is 

# = 4 = 0.88. (6) 

^ s.o. ^'-'L 

This ratio is very well reproduced in the corresponding spectra of Figs. 2 and 
3. Concerning the ratios of the spin-orbit splitting of the proton and neutron 
to the S+ and S° hyperon, respectively, a simple expression similar to the 
one of Eq. cannot be given because of the different couplings. However, 
in the corresponding spectra of Figs. 2 and 3 we found a value of about 
~ 0.34. 

In Table III we compare various S^, and T^^ levels for the |?Si, f Ca, 
|PZr, and I^^Pb hypernuclei. Of course, the Coulomb force plays an impor- 
tant role: comparing S~ with TP levels, the atomic states disappear, whereas 
for the S"*" states the baryonic potential has to overcome the Coulomb repul- 
sion with the effect that only the deep lying states survive. For symmetric 
(A^ = Z) hypernuclei, where the p-meson field is weak (it is nonzero because 
the proton and neutron density distributions differ due to the Coulomb in- 
teraction), we found the Coulomb shifts between E~ and or TP and S"*" 
states almost identical to the corresponding neutron-proton shifts in "nor- 
mal" symmetric nuclei. The situation is somewhat different for hypernuclei 
with a neutron-excess because then the effect of the p-meson field is not 
negligible anymore: For the T" (S+) the p-meson adds (subtracts) from the 
isoscalar part of the timelike repulsive vector field. The S° does not couple 
to the p field at all. To get an idea of the impact of the p-meson field we 
recalculated the asymmetric hypernuclei |PZr and |P^Pb with the p-T cou- 
pling switched off. For |PZr we found the S^(S+) states stronger (weaker) 
bound by about 2.1-3.5 MeV; the corresponding range for |°^Pb is 4.4-6.9 
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MeV. The fact that such relatively large ranges occur can be understood in 
terms of the rms-radii: The lower bounds are for weakly bound states with 
large rms-radii (e.g. the 1(79/2 S~ state with r^ms = 5.40 fm for |PZr and 
the 3pi/2 S~ state with Vrms = 6.97 fm for |°^Pb), while the upper bounds 
correspond to deep lying states with small rms-radii (e.g. the lsi/2 S~ states 
with Trms = 3.09 fm for |?Zr and Vrms = 3.68 fm for |°^Pb; the values for the 
rms-radii are from the calculations with the p field switched on). As one can 
see, the influence of the p-meson field, whose range is determined by its mass 
rrip, weakens with increasing radial distances. (For neutrons and protons the 
situation is similar to the S~ — S+ pair but the effective p-nucleon coupling 
is weaker.) 

Hence, the quantum hadrodynamical treatment of hypernuclei offers, by 
the possible inclusion of the p-meson field, a natural way to incorporate an 
isospin dependence into the S potential (i.e. a Lane potential), which was 
pointed out by Dover in Ref. ^ to be one of the most important questions 
of hypernuclear physics. 

Finally, we show in Figs. 4-6 the single-particle energies of the 
and S+ hyperons, respectively, versus A~'^^^, with, A the mass number of the 
nuclei. For the S~ (Fig.4) the attractive Coulomb potential alone, irrespec- 
tive of the strength of the short-range nuclear potential, is enough to bind 
S~ states. Some of these states (the most bound) are such that the rms-radii 
of the S~ wavefunctions are essentially inside the nucleus: Looking at |P_Ca 
we found for the rms-radii of the plotted ss-, ps-, ds-, and /s- states the 
values of Vrms = 2.65, 3.31, 3.89, and 4.76 fm, respectively. Therefore, com- 
pared with the experimantal charge rms-radius of ~ 3.48 fm for ^°Ca, the 
(is- and /s- states should be called atomic ones, while for the other levels 
an identification as hypernuclear states seems to be more appropriate. 

For the charge neutral S° (Fig. 5) there is no Coulomb attraction and 
therefore the number of bound states decreases. The value of -28 MeV for 
y4~^/^ = 0.0 represents the binding energy of the lowest E° level in nuclear 
matter, which is under the assumption of an universal hyperon coupling the 
same as for the A particle |]T], ^. As expected, the pattern of states shows 
the standard behaviour as for A hypernuclei 0. 

Turning finally to the discussion of the A-dependence of the S+ levels 
(Fig. 6) the situation becomes more complicated. Of course, now the Coulomb 
force is repulsive and the number of bound states further decreases compared 
with S~ and S° hypernuclei. But to get a full understanding of Fig. 6, it 
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seems necessary to consider the various contributions to the nuclear potential 
in view of their range, which is determined by the corresponding meson 
mass, respectively. For example, the fact that the binding of the state 
increases going from A = 90 to A = 208, while there is an opposite effect 
for the Ss+ and ps+ states, may be attributed to the larger 0?^+ rms-radius 
(e.g. 5.63 fm instead of 4.55 fm and 5.17 fm for the ss+ and ps+ levels in 
|°+Pb, respectively). The wavefunction is located at large radial distance, 
where the impact of the attractive cr-meson (m^ = 500 MeV) increases locally 
compared with the repulsive cu-meson (m^ = 783 MeV) due to their different 
ranges. Seemingly, the various states show a behaviour similar to the one 
found within nonrelativistic Hartree-Fock Skyrme calculations for protons in 
"normal" nuclei, where the binding of the d levels strongly increases, whereas 



the s states nearly stay constant when going from A =90 to A =208 [T^ . 

In the present calculations the broadening of the S hyperon states due to 
their decay to the A was neglected. In principle, the model can be extended 
to include the decay by introducing appropriate EA-vertices. An important 
aspect of such an extended approach would be the possibility of investigating 
the decay of strange particles in the nuclear medium. In order to estimate 
the effects due to the conversion SA^ AN the results of nonrelativistic 
potential models [0 may be taken as a guideline. In such approaches the 
decay is described schematically by an absorption potential for hyperons. The 
imaginary part of the self-energy effectively lowers the binding energy which 
can be understood in terms of the pole structure of the baryon propagator. 
Qualitatively a similar effect has to be expected also in a covariant description 
including the decay of the E hyperons. Thus the present results are likely to 
give lower bounds for the binding properties of the E particles in hypernuclei. 

In conclusion, we performed relativistic mean-field calculations of A and 
E hypernuclei using an interaction that considers neutron-star masses, the A 
binding in saturated nuclear matter, and experimental A single-particle lev- 
els. Concerning the E couplings we assumed an universal hyperon coupling; 
i.e. all hyperons in the lowest octet couple to the meson fields as the A. 
We employed the so-called expectation-value method whose reliability was 
found to be sufficient compared with fully self-consistent RHA calculations. 
Analyzing the hypernuclear spectra, we devoted special attention to the spin- 
orbit interaction for hyperons and, in the case of E hypernuclei, to the isospin 
dependence of the interaction. These two features are of particular interest 
in the current discussion concerning hyperon potentials in nuclei and are nat- 
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urally incorporated into the relativistic quantum hadrodynamical model we 
used. In addition we found by comparison with corresponding nonrelativistic 
results that for strongly bound states the impact of the T,N — > AN decay 
width on the bindings seems to be negligible. 

In the future it would be very valuable from the point of view of dense 
matter properties, and especially the structure of neutron-stars, to have the 
assumption of an universal hyperon coupling confirmed by detailed precision 
experiments on E hypernuclei, and we hope that these calculations may 
possibly be of assistence as well as a stimulus to such experiments and the 
development of the necessary facilities. 
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Table captions 



table I: Parameters of the two forces considered in the text. In both cases 
the nucleon mass M = 938 MeV. For saturated nuclear matter the set 
11 1^ yields energy per particle E/A = —16.3 MeV, density po = 0.153 
fm~^, incompressibility K = 300 MeV, effective mass M* /M = 0.7 and 
symmetry energy coefficient agym = 32.5 MeV. 
Cf = gf {M/rriif ,Xi = gm/gui = (T,uj,p;H = A, E. 

table II: Comparison of our H* results for various A levels in ^Ca. and ^°^Pb 
with the fully self-consistent relativistic Hartree results from Ref. [Q]. 
All quantities are in MeV. 

table III: Various S^, and S+ levels in f Si, f Ca, |°Zr, and |°spb. All 
quantities are in MeV. 
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Table I 





m 


Hi 


(MeV) 


499.31 


500 


(MeV) 


780 


783 


mp (MeV) 


763 


770 


Ma (MeV) 


1116.08 


1115 


Ms (MeV) 




1190 




348.26 


266.40 


CI 


229.29 


161.53 


CI 


148.92 


99.67 


b X 10^ 


2.2847 


2.947 


c X 10^ 


-2.9151 


-1.070 




0.464 


0.600 




0.481 


0.653 


Xp 




0.600 
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Table II 




1(^3/2 


-2.63 


-1 


lc^5/2 


-3.76 


-2 


lPl/2 


-10.93 


-8 


lP3/2 


-11.61 


-9 


lSl/2 


-19.43 


-16 



17 -15.78 -15.03 

08 -16.12 -15.32 

75 -20.38 -19.42 

38 -20.51 -19.54 

90 -24.19 -23.23 
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Table III 

















Level 


E- 




E+ 


E- 


EO 


E+ 


1/5/2 


- 


- 


- 


-2.31 


- 


- 


1/7/2 


- 


- 


- 


-3.43 


- 


- 


2si/2 


-5.19 


-0.43 


- 


-9.73 


-2.57 


- 


la!3/2 


-5.10 


- 


- 


-10.29 


-3.00 


- 


1^5/2 


-6.19 


-0.87 


- 


-11.31 


-3.99 


- 


lPl/2 


-13.97 


-7.96 


-2.07 


-18.79 


-10.82 


-3.00 


lP3/2 


-14.71 


-8.69 


-2.79 


-19.37 


-11.42 


-3.60 


lSl/2 


-23.22 


-16.63 


-10.11 


-27.10 


-18.54 


-10.07 






|PZr 






208pb 




1/5/2 


-13.35 


-4.03 




-26.47 


-12.40 




1/7/2 


-14.30 


-5.00 




-26.87 


-12.87 




2Si/2 


-19.35 


-9.28 




-30.79 


-15.76 


-0.73 


14/2 


-20.29 


-10.52 


-0.87 


-31.57 


-16.99 


-2.70 


14/2 


-20.88 


-11.15 


-1.53 


-31.76 


-17.24 


-3.02 


lPl/2 


-26.90 


-16.61 


-6.47 


-36.25 


-21.06 


-6.22 


lP3/2 


-27.16 


-16.91 


-6.80 


-36.31 


-21.15 


-6.36 


lSl/2 


-32.94 


-22.05 


-11.31 


-40.44 


-24.48 


-8.94 
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Figure captions 



figure 1: Hyperon self-energy contributions gHa(^ (dotted lines), Qhw'^^ 
(dashed lines), (y'spPo (long-dashed lines), and eA^ (dot-dashed lines) 
for the nuclei ^^Si, ^°Ca, ^°Zr, and ^''^Pb. The nonrelativistic reduction 
of the hyperon potential is (for A and TP entirely, for mainly) given 
by the difference Qhu^'^ — QHaO', which is represented by the solid lines, 
respectively, if = A, S in the hyperon coupling constants. 

figure 2: The calculated proton, S"*", neutron. A, and single-particle 
spectrum for ^'^Ca (parameter set II of Table I). 

figure 3: Same as Fig.2 for ^°^Pb. 

figure 4: Single-particle energies versus A'"^^^ for S^. For each angular 
momentum the lowest lying state is plotted, respectively. The dashed 
lines are added to guide the eye. 

figure 5: Same as Fig.4 for S°. The value of -28 MeV for A^^/s = o.O 
represents the binding energy of the lowest S° level in saturated nuclear 
matter under the assumption of an universal hyperon coupling [Q, |]. 

figure 6: Same as Fig.4 for S^. 
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